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Abstract. The y-process in core-collapse and/or type la supernova explosions is 
thought to explain the origin of the majority of the so-called p nuclei (the 35 proton-rich 
isotopes between Se and Hg). Reaction rates for 7 -process reaction network studies 
have to be predicted using Hauser-Feshbach statistical model calculations. Recent 
investigations have shown problems in the prediction of a-widths at astrophysical 
energies which are an essential input for the statistical model. It has an impact on the 
reliability of abundance predictions in the upper mass range of the p nuclei. With 
the measurement of the 164 , 166 Er(a,n) 167,169 Yb reaction cross sections at energies 
close to the astrophysically relevant energy range we tested the recently suggested 
low energy modification of the a+nucleus optical potential in a mass region where 
7 -process calculations exhibit an underproduction of the p nuclei. Using the same 
optical potential for the a-width which was derived from combined 162 Er(a,n) and 
162 Er(a,y) measurement makes it plausible that a low-energy modification of the 
optical a+nucleus potential is needed. 
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1. Introduction 

Contrary to the synthesis of the light nuclei - where charged particles are playing a 
particularly important role [ 1 ] - two different neutron-capture processes, the so-called s 
[2] and r [3] processes, are required to produce the bulk of naturally occurring nuclides 
above Fe. These two processes were found to be unable, however, to create 35 neutron- 
deficient, stable, rare isotopes between ' 4 Se and 19 b Hg, which were called u p nuclei” or 
“excluded isotopes” [4, 5]. Photodisintegration of stable nuclei either in the outer shells 
of massive stars during a core-collapse supernova explosion [ 6 , 7] or during the explosion 
of a White Dwarf [ 8 ] has been suggested as a production mechanism for most of these 
nuclei ||. Such a so-called 7 -process commences at 2-3 GK temperatures by sequences 
of (y,n) reactions, which are replaced with (y,p) and ( 7 , a) reactions when reaching 
sufficiently neutron deficient nuclides in an isotopic chain. The whole process lasts only 
a few seconds, before the environment becomes too cool for the photodisintegrations. 
Recent calculations, however, are under-predicting the heavy p-nucleus abundances in 
the mass region between 150 < A < 165 [10, 11, 12 ], It remains unclear whether this 
deficiency is due to nuclear cross sections, stellar physics, or if alternative or additional 
processes have to be invoked [ 11 ]. 

In order to resolve this ambiguity, on one hand, improvements on the description 
of the astrophysical conditions under which the process takes place (seed isotope 
abundances, peak temperatures, time scale, etc.) are needed. On the other hand, 
uncertainties are introduced into the calculations by nuclear physics input, in particular 
by the reaction rates. The 7 -process models require the use of huge reaction networks 
including ten thousands of nuclear reactions and the rates of these reactions at a given 
stellar temperature have to be known. The reaction rates are generally taken from 
calculations using the Hauser-Feshbach (H-F) statistical model [13, 14, 15, 16, 17]. The 
accuracy of the H-F predictions mainly depends on the adopted nuclear models for the 
proton-, neutron-, a-, and 7 -widths. The predictions can be tested for only a limited 
number of reactions, since experimental data in the relevant mass and energy range are 
very scarce [18]. Although photodisintegrations play a role in the 7 -process, their rates 
are computed from capture rates using detailed balance. Experimental information on 
quantities (such as particle- and 7 -widths) for the calculation of the stellar rates (which 
cannot be directly measured, with a few exceptions[19]) can be obtained from the study 
of capture reactions. This approach is not only technically less challenging, but also 
provides more relevant astrophysical information than the direct study of the 7 induced 
reactions [17, 20, 21, 22, 23]. 

Reactions with a particle absorbtion or emission in the 7 -process are mainly 
relevant above the neutron shell closure at N = 82 [11, 23]. Previous 7 -process 
related experiments found discrepancies between measured and predicted a-induced 

|| The origin of few isotopes — such as the 164 Er, 152 Gd and 180 Ta nuclei — is still under debate. 
Recently it was shown that there is a large s process contribution to their observed abundances [9, 11]. 
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Figure 1 . (Color online) XRF (A) and PIXE (B) spectra used to determine the 
impurities and the thickness of the erbium targets. The peaks used for the analysis 
are marked. Peaks belonging to the Rh anode of the X-ray tube, to impurities in the 
target and/or the backing are indicated, too. For details, see text. 


cross sections, which suggested low-energy modifications of the o+nucleus optical 
potential [24, 25, 26, 27, 28, 29] or the need to modify the reaction model [30]. An 
important recent break-through was the simultaneous, high-precision measurement of 
the 162 Er(a, 7 ) and 162 Er(ct,n) reaction cross sections in a wide energy range below the 
Coulomb barrier [29]. It was crucial for the theoretical interpretation that (a, 7 ) data 
were taken even below the (a,n) threshold. Thus, for the first time it was possible to 
show that the cc+nucleus optical potential requires an energy-dependent modification 
consistently and unambiguously within the same measurement at high masses. The aim 
of the present work is to provide further data in the same energy region for other Er 
isotopes and to test whether they can be described by the same approach as the 162 Er 
data. 

The paper is organized as follows: Experimental details are described in Sec. 2; the 
resulting cross sections are presented in Sec. 3, where also a comparison with theoretical 
predictions is made; Sec. 4 provides the conclusions and a short summary. 

2. Experimental technique 

The experiment was carried out at the Institute for Nuclear Research, Hungarian 
Academy of Sciences (MTA Atomki, Debrecen, Hungary). The element Er has six 
stable isotopes with mass numbers A = 162, 164, 166, 167, 168 and 170. Alpha induced 
cross sections on the p nucleus 162 Er were already measured [29], the aim of the present 
work was to study the 164 Er(a,n) 16 'Yb and 166 Er(a,n) 169 Yb reactions. To determine the 
a-induced reaction cross sections the activation method was used; the 7 -rays following 
the electron capture decay of the Yb reaction products were detected using a Low Energy 
Photon Spectrometer (LEPS) and a coaxial HPGe detector. The following paragraphs 
provide a detailed description of the experimental setup. 
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Figure 2. Particle spectrum used to monitor the target thickness (A). The peaks 
corresponding to alpha backscattering on erbium, on the aluminum backing and on 
target impurities are indicated. On the right side the backscattering yield as the 
function of the current integrator counts is shown (B), a linear function was fitted to 
the data to prove that the target is stable. For more details see the text. 


2.1. Target production, characterization 

The targets were made by vacuum evaporation of Er 2 03 powder with natural isotopic 
composition ((1.61 ± 0.03) % 164 Er and (33.6 ±3.5)% 166 Er, respectively) and enriched 
to (8.1 ±0.2)% in 164 Er ( 166 Er content: (33.10 ± 3.4) %) onto 2 /irn thick, high purity 
A1 foils. To check the enrichment provided by the supplier, at E Q = 15.5, 16.5 and 
17.0 MeV the irradiation was carried out using both natural and enriched targets and 
the derived cross sections were found to be consistent. 

The Er targets were produced using reductive evaporation technique, the Er 2 03 
powder was mixed with Zr powder and placed into a carbon crucible heated by electron 
beam. During the evaporation, the distance between the A1 backing and the evaporation 
boat was 8 cm, using this distance, uniform targets could be produced. Care was taken to 
limit the target’s Zn contamination — identified as the main contaminant of the backing 
and the Zr powder used for the reduction — since during the a irradiation 67 Ge could be 
produced via the 64 Zn(a,n) 67 Ge reaction. This isotope (and its daughter) was identified 
as the main disturbing activity since during the /? decay a 208.95 keV 7 ray is emitted 
which is only 1.15 keV above the 207.80 keV 7 of our interest (see Table 1). As a first 
step, the quality of the targets were tested by X-ray fluorescence spectroscopy (XRF) 
[31]. The targets were excited by hltered primary radiation (with a 2 mm A1 filter) of 
an Oxford Instruments X-ray tube (XTF 6000 BR) having Rh target. It was operated 
at 30 kV voltage and 700 pA electron beam current. The measurement time was 1000 s 
in each case. The characteristic X-rays were detected in atmospheric pressure condition 
enabling the detection of trace elements down to potassium with high sensitivity. The 
evaluation of the spectra was carried out with the AXIL software [32], 

The absolute thickness of the targets were measured by weighing and using the 
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Figure 3. Off-line 7 spectra, taken after the irradiation of an erbium target with 
14.5 MeV a beam. The upper figures show the 7 lines used to determine the cross 
section of the 164 Er(a,n) 16 'Yb reaction: in figure A the 7 transitions emitted during 
the decay of the 16 'Yb are indicated. During the decay of the 167 Tm nucleus — the 
daughter of the 167 Yb reaction product — a 207.80 keV 7 ray is emitted, figure B shows 
this transition (together with the 208.95 keV background transition emitted during 
the decay of the 67 Ge produced via 64 Zn(a,n) 67 Ge reaction) and figure C presents its 
decay curve. The 7 -lines used to determine the cross section of the 166 Er(a,n) 169 Yb 
are marked on the lower figure. For further details, see text. 


PIXE (particle induced X-ray emission) technique [33]. The PIXE chamber is installed 
on the left 45° beam line of the 5 MV VdG accelerator of Atomki. A more detailed 
description of the PIXE setup can be found in Ref. [34], A homogeneous beam of 
2 MeV protons with a 5 mm diameter and a 4-5 11 A current was used for the thickness 
measurement. The total collected charge in the case of each target was about 2 fiC. The 
measured spectra were fitted using the PIXEKLM code [35]. Typical XRF and PIXE 
spectra can be seen in Figure 1 where peaks used for the analysis are marked. From 
the analysis of the PIXE spectra the erbium target thickness was determined and found 
to be in agreement with the weighing within maximum of 4.4 %. This 4.4 % — as a 
conservative estimate — has been adopted as the uncertainty of the target thickness 
determination. 

The thickness of the targets varied between 114 and 238/ug/cm 2 , corresponding to 
an a energy loss of AE 15 keV (at E a = 17MeV) and AE fa 39 keV (at E Q = 12 MeV), 
respectively. The a energy loss was calculated using the SRIM [36] code. Thicker targets 
were used at low bombarding energy, where the cross section and the corresponding 7 
yield is smaller. 
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Table 1. Decay parameters of the 164 Er(a,n) 167 Yb and 166 Er(a,n) 169 Yb reaction 
products taken from the literature [37, 38]. 


Reaction 

Threshold Product 
(MeV) nucleus 

Half-life 

(hour) 

7 -ray 

energy (keV) 

Relative 7 -intensity 
per decay (%) 

164 Er(a,n) 167 Yb 

11.27 

167 Yb 0.292 ±0.003 

106.2 

22.44 ±1.35 





113.3 

55.08 ±2.97 





176.2 

20.40 ±0.80 



167 Tm 

222 ±4.8 

207.8 

41.54 ±8.00 

166 Er (a,n ) 169 Yb 

10.44 

169 Yb 

768.4 ±1.2 

63.1 

43.62 ±0.23 





109.8 

17.39 ±0.09 





130.5 

11.38 ±0.05 





177.2 

22.28 ± 0.11 





197.9 

35.93 ±0.12 


2.2. Irradiation and 7 counting 

The erbium targets were irradiated with a beams from the MGC 20 cyclotron of Atornki. 
Energies were in the range of 12<E a <17MeV, covered in steps of about 0.5 MeV. 
After the beam-defining aperture, the chamber was insulated and a secondary electron 
suppression voltage of -300 V was applied at the entrance of the chamber. From the 
last beam-defining aperture the whole chamber served as a Faraday cup. The collected 
charge was measured with a current integrator, the counts were recorded in multichannel 
scaling mode, stepping the channel in every minute to take into account the possible 
changes in the beam current. The length of the irradiation were between 7 and 25 hours 
corresponding to number of incident a particles in each irradiation between 0.9 x 10 1 ' 
and 6.1 x 10 17 . 

Before the experiment several beam tests were performed to determine the 
maximum tolerable a beam current. These tests showed that there was no deterioration 
of the targets using an a beam current less than 2.5 /iA. The target thickness was 
monitored by measuring the yield of the back-scattered a particles, for this purpose an 
ion implanted silicon detector built into the chamber at 9 = 165° relative to the beam 
direction was used. The spectra of this detector were recorded regularly and the ratio of 
the alpha backscattering yields divided by the corresponding current integrator counts 
were fitted using a linear function. It was found that the target-losses are typically in 
the order of 0.6% and the maximum target loss was 1.2 ± 0.15%. Figure 2 shows a 
typical alpha backscattering spectrum (A) and the corresponding stability curve (B) 
with the fitted linear function. 

After the irradiations, 0.25 h waiting time was used in order to let short-lived 
activities, which would impact the quality of the measurement, decay. The duration of 
the 7 countings were about 200 - 300 h in case of each irradiation. The decay parameters 
of the reaction products are given in Table 1. To determine the 164 Er(a,n) 167 Yb reaction 
cross section first the decay of the produced 16 ' Yb (T 72 = 0.292 ± 0.003 h) was followed, 
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Table 2. Measured cross sections of the 164 Er(a,n) 167 Yb reaction. 


.Eiab E c m . Cross section 

[MeV] [MeV][mbarn] 


13.5 

13.17 

± 

0.04 

0.069 ± 

0.006 

14.0 

13.66 

± 

0.04 

0.208 ± 

0.018 

14.5 

14.14 

± 

0.07 

0.665 ± 

0.058 

14.51 

14.16 

± 

0.15 

0.779 ± 

0.079 

15.0 

14.63 

± 

0.08 

1.55 ± 

0.15 

15.5 

15.12 

± 

0.05 

4.62 ± 

0.41 

16.0 

15.62 

± 

0.08 

10.3 ± 

0.92 

16.5 

16.09 

± 

0.08 

18.6 ± 

: 1.6 

16.52 

16.11 

± 

0.05 

20 ± 

1.8 

17.0 

16.59 

± 

0.05 

29 ± 

2.7 


the yield of the 106.2 keV and 113.3 keV 7 transitions were measured. Moreover, the 
167 Tm nucleus (Ti / 2 = 222 ± 4.8 h), the daughter of the produced unstable 167 Yb, 
decays by electron capture to 16 'Er with emission of 207.8keV 7 ray, which was also 
used to determine the reaction cross section. The determination of the 166 Er(a,n) 169 Yb 
reaction cross section is based on the measurement of the 7 rays emitted during the (3 
decay of the produced 169 Yb (e.g. see figure 3). 

For the 7 counting a low-energy photon spectrometer (LEPS) was used. In order 
to reduce the laboratory background a multilayer quasi 47 T shield was built around the 
LEPS detector including an inner 4 mm thick layer of copper, a 2 mm thick layer of 
cadmium, and an 8 cm thick outer lead shield [40]. To limit the systematic uncertainty 
the activity of several samples irradiated at 15.0, 16.0 and 17.0 MeV were measured 
using a coaxial HPGe detector with 100 % relative efficiency, placed in a low background 
shielding, too. The resulted cross sections are within 5.2% compared to the ones based 
on the counting with the LEPS detector. 

I 11 the case of both detectors, the distance between the target and the end cup was 
1 cm. The detector efficiencies for both detectors used in the present work had to be 
known in this close geometry with high precision. A similar procedure was used to derive 
the photopeak efficiencies in the case of the coaxial HPGe and LEPS detectors. First the 
absolute detector efficiency was measured in far geometry: at 15 cm distance from the 
surface of the detector, using calibrated 57 Co, 133 Ba, 13 'Cs, 152 Eu, and 241 Am sources. 
Since the calibration sources (especially 133 Ba and 152 Eu) emit multiple 7 radiations from 
cascade transitions, in the close geometry a strong true coincidence summing effect is 
expected, resulting in an increased uncertainty of the measured efficiency. Therefore, no 
direct efficiency measurement in close geometry was carried out. Instead, the activity 
of several irradiated erbium targets were measured in both close and far geometries. 
Taking into account the time elapsed between the two countings, a conversion factor of 
the efficiencies between the two geometries could be determined and used henceforward 




Table 3. Measured cross sections of the 166 Er(a,n) 169 Yb reaction. 


.Eiab E c , m . Cross section 

[MeV] [MeV][mbarn] 


12.00 

11.70 

± 

0.07 

0.0028 

± 

0.0004 

12.50 

12.18 

± 

0.07 

0.011 

± 

0.001 

13.00 

12.68 

± 

0.04 

0.027 

± 

0.003 

13.50 

13.17 

± 

0.04 

0.104 

± 

0.009 

14.00 

13.66 

± 

0.04 

0.264 

± 

0.023 

14.51 

14.16 

± 

0.15 

0.89 

± 

0.08 

15.00 

14.63 

± 

0.08 

2.26 

± 

0.19 

15.01 

14.64 

± 

0.08 

2.06 

± 

0.17 

15.50 

15.12 

± 

0.05 

5.21 

± 

0.44 

16.00 

15.62 

± 

0.08 

11.3 

± 

0.96 

16.50 

16.11 

± 

0.05 

21.2 

± 

1.79 

17.00 

16.59 

± 

0.05 

37.2 

± 

3.17 

17.01 

16.61 

± 

0.05 

35.9 

± 

3.17 


in the analysis. 

3. Results and discussion 

The measured (a,n) cross section values are listed in Tables 2 and 3. The quoted 
uncertainty in the E c _ m . values corresponds to the energy stability of the a beam and to 
the uncertainty of the energy loss in the target, which was calculated using the SRIM 
code [36]. The uncertainty of the cross sections is the quadratic sum of the following 
partial errors: efficiency of the detectors ( 6 % (HPGe) and 5% (LEPS), respectively), 
number of target atoms (4.4%), current measurement (3%), uncertainty of decay 
parameters (<8.0%), and counting statistics (0.4% - 7.4%). Some irradiations were 
repeated at the same energies. The cross sections were then derived from the averaged 
results of the irradiations weighted by the statistical uncertainty of the measured values. 

The prediction of averaged reaction channel widths is crucial for the Hauser- 
Feshbach model of nuclear reactions [13, 17]. In this well-established model it is assumed 
that a large number of unresolved resonances is located around the formation energy 
Ec = E c . m . + F“ p of the compound nucleus (in our case 168 , 170 Yb), with F“ p being 
the a-separation energy in the compound. Instead of individual resonance widths, an 
average over all resonance widths is used, leading to averaged neutron-, proton-, a-, 
and 7 -widths (further widths are negligible at the investigated energies). These widths, 
describing the formation and decay of the compound nucleus states are predicted from 
further models. Following [22, 23], the sensitivity 

V da 
8 ~ a dT 


( 1 ) 
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Figure 4. (Color online) Sensitivities of the investigated reaction cross sections (A: 
164 Er(a,n) 167 Yb, B: 166 Er(a,n) 169 Yb) to the a, 7 and neutron widths. For details, see 
text. 


describes the change in the cross section da when a width T is changed by dr. 
Not all the widths are important at each energy. The sensitivity of the calculated 
cross sections to variations of the different widths is plotted in Figure 4 (the proton 
width is not shown as the cross section is insensitive to it across all shown energies). A 
sensitivity of unity means that the cross section is changed by the same factor / as the 
input parameter, whereas the cross section does not change when the sensitivity is zero. 
As can be seen the cross sections depend dominantly on the a-width at the upper end of 
the measured energy range. This is because the averaged a-width remains much smaller 
than the neutron width at energies well above the reaction threshold. Towards lower 
energy approaching the reaction threshold, the neutron width strongly decreases and 
the cross sections become increasingly sensitive to the neutron- and 7 -widths in addition 
to the a-widths. It was already pointed out in [28] that this behavior does not allow 
for an unambiguous identification of the source of discrepancy closer to the threshold 
from (ct,n) data alone, if such a discrepancy between cross section measurements and 
prediction were found. 

A similar situation as in [28] is encountered here. Figure 5 shows a comparison 
of our data to Hauser-Feshbach predictions performed with the code SMARAGD [41], 
In the case of the 166 Er(a,n) 16 Wb reaction, figure 5 also includes the previous data by 
[28] which are in good agreement with the data obtained in this work. While there 
is excellent agreement of the data with the standard calculation (marked by the full 
line and denoted by “McF-S” in the figures) at higher energy, an increasing deviation 
is found towards lower energy for both reactions. This indicates that the a-width is 
described well at the higher energies. The main ingredient for the calculation of the 
averaged a-width is the a+nucleus optical potential, for which the global potential by 
McFadden and Satchler [42] was used. The discrepancies towards lower energy cannot 
be unambiguously assigned to a misprediction of a certain width when solely using the 
present (a,n) data. 
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Figure 5. (Color online) Measured reaction cross sections of the 164 Er(a,n) 167 Yb 
(A) and 166 Er(a,n) 169 Yb (B) reactions compared to Hauser-Feshbach predictions (see 
text). 


A recent measurement provided cross sections for both 162 Er(a,n) 165 Yb and 
162 Er(a, 7 ) 165 Yb reactions, with (a, 7 ) data also below the (a, 11 ) threshold [29]. Due 
to the additional (a, 7 ) data it was possible to show in [29] that an energy-dependent 
modification of the optical a+nucleus potential is required at low energies. Additionally, 
close to the (a,n) threshold the ratio of the neutron- to 7 -width had to be increased 
by about 40% (without the possibility to assign this change to either width). For 
comparison, in Figure 5 we show also cross sections obtained from calculations using 
the same modifications of widths as in [29] (dashed lines labelled “mod”). Again, good 
agreement is obtained without further changes. It should be noted that a modification 
of the neutron- to 7 -width ratio only affects the cross sections very close above the 
reaction thresholds (up to about 0.7 MeV above threshold) because of the strongly 
decreasing sensitivity to those widths with increasing energy, which can be seen in 
Figure 4. Therefore, the changes in the cross sections in the energy ranges covered by 
the present measurements, and shown in Figure 5, are solely due to the modification of 
the a widths. 

4. Summary 

For the first time, (a,n) cross sections of the p nucleus 164 Er were measured at energies 
well below the Coulomb barrier, using the activation technique. Furthermore, the 
166 Er(o;,n) 169 Yb reaction was studied at energies overlapping with the results of a recent 
measurement performed by [28]. The agreement between the two data sets is almost 
perfect, this fact strengthens further the conclusions of [28]. 

The experimental data were compared to statistical model predictions. Good 
agreement was found at the upper end of the measured energy range when using the 
global a+nucleus optical potential by [42] but deviations between data and calculations 
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were found towards the reaction threshold in both reactions. Using the same set of 
parameters as determined in a previous, simultaneous measurement of 162 Er(a, 7 ) and 
162 Er(a,n) [29] works well also in the present cases. This confirms the necessity of an 
energy-dependent modification of the a+nucleus optical potential at very low energies. 

The present measurements enlarge the database for deriving an improved global 
a+nucleus optical potential at very low energies by future theoretical models. The 
relevant astrophysical energy range, however, is below the reaction threshold and 
therefore further information on the energy dependence of the a widths at astrophysical 
energies requires additional experiments. 
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